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Current- voltage curves of small area hysteretic intrinsic Josephson junctions, for which the Joseph- 
son energy Ej — hj c /2e is of order of thermal energy kT, are investigated. A non-monotonic 
temperature dependence of the switching current is observed and explained by thermal phase fluc- 
tuations. At low temperatures premature switching from the superconducting into the resistive state 
is the most important effect of fluctuations. At high temperatures only a single resistive branch is 
observed. At the cross-over temperature a hysteretic phase-diffusion branch exists. It shows the 
importance of a frequency-dependent impedance of an external circuit formed by the leads. 

PACS numbers: 74.40.+k, 74.50.+r 



-a 
c 

o 
o 



> 

(N 

O 
m 
o 

-t— > 



i 

-a 
c 

o 
o 



I. INTRODUCTION 

In the strongly anisotropic cuprate supercon- 
ductors, such as Bi 2 Sr 2 CaCu 2 0-8+,5 (BSCCO) or 
Tl 2 Ba 2 CaCu 2 8+5 (TBCCO), the Cu0 2 layers together 
with the intermediate material form a stack of Josephson 
junctions. In the presence of a bias current perpendicular 
to the layers each junction of the stack is either in the 
resistive or in the superconducting state leading to the 
well- know n multibranch structure of the IV-curves, see 
Ref. lHI313 and references therein. 

The effect of thermal noise on the properties of Joseph- 
son junctions made from low-temperature superconduc- 
tors, and in particular on the critical current and voltage- 
current curves, has been investigated experimentally and 
theoretically in many publications, see Ref. 5 6 7 8. For 
high-temperature superconductors considerable interest 
has been attracted to the problems of therma l fluctu- 
ations in the vortex state, see Ref. lPllflllllUlll but, 
to our knowledge, no detailed study of the influence of 
thermal fluctuations on the critical current and current- 
voltage curves of mesoscopic intrinsic Josephson junc- 
tions has been carried out. Some fluctuation effects are 
reported recently in Ref.ll4l 

In our paper we present the results of experimen- 
tal and theoretical investigations of small area intrin- 
sic Josephson junctions, when the Josephson energy 
Ej{T) = J c (T)/2e is of the order of kT (k is the Boltz- 
mann constant). For temperatures T > T* thermal 
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Figure 1: Estimation of crossover temperature T* from the 
condition E.j(T*) m kT* for mesas of different area. 



fluctuations become important. This cross-over tem- 
perature can be estimated with help of Fig.^ where 
the Ambegaokar-Baratoff temperature dependence of the 
critical current density j c {T) is plotted together with the 
linear functions 2ekT/hA (marked as kT = Ej in the fig- 
ure) for mesas with different area A. T* can be read-off 
from the crossing points of the curves. 

We discuss in the following junctions with cross-over 
temperatures in the range of 10-60 K (A~ 0.5x0.5/xm). 
Note that for all our samples the charging energy 
E c = e 2 /2C is still small compared to Ej, such that 
quantum effects (Coulomb blockade and quantum tun- 
neling) are not essential. From the behaviour of the 
current- voltage characteristics of the junctions three tem- 
perature regimes can be distinguished. At low temper- 
atures, T < T* , straight lines with hysteretic jumps 
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to other resistive branches are observed (see Fig.0 and 
Fig.|8J|. At temperature T ~ T* the current voltage 
curves are bended near the critical current (like in Fig.0] 
and Fig. lllf) . but still some hysteresis is found. At high 
temperatures, T > T* , no more jumps in the current 
voltage curve occur and no real critical current exists, 
but still a maximum in the curvature can be distinguished 
(see Fig.|SJl, which can be regarded as maximum of the su- 
percurrent. In order to describe all temperature regimes 
together we introduce the concept of a " switching" cur- 
rent, which will be defined below. It replaces the critical 
current in the non-hysteretic high-temperature regime. A 
plot of the switching current as function of temperature 
shows a rather non-monotonic behaviour, which is very 
different from the AB-curve found for the critical current 
of large junctions. We will show below that the temper- 
ature dependence of the switching current, in particular 
the sharp drop, can be explained by phase diffusion pro- 
cesses. This will be done with help of an extended RSJ 
model containing thermal current noise and the coupling 
to an external RC circuit (Fig.[2J). The latter is necessary 
to obtain the observed hysteretic behaviour for the phase 
diffusion branch at intermediate temperatures T ~ T* . 
In our case the external RC circuit is formed most prob- 
ably by a large capacitance between the ground super- 
conducting layer and the upper electrode together with 
the resistance of the contact between mesa and upper 
electrode. 

The paper is organized as follows. In section [H] we 
discuss typical voltage-current curves at different tem- 
peratures and give a brief summary of the theory used 
in calculations. In section IIIII sample preparation and 
experimental set-up are described. The main section ITVl 
is devoted to experimental results and their discussion. 
Here also a discussion of magnetic field effects is included. 




Figure 2: RSJ model with additional frequency dependent 
damping modeled by an external RC circuit. 



{j T (T)jT(r')) = 2 7 6(T-r'), 



(4) 



where r = uj c t, j = J j J Cl and lu c — 2eRJ c /h is the 
characteristic frequency. 

Here we introduce two main parameters: 
(i) the McCumber parameter 



(5) 



where u> 2 — 2eJ c /hC is the Josephson plasma frequency; 
(ii) the fluctuation parameter (normalized temperature) 



kT 2ekT 



TlJ c 
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Extended RSJ model 



II. OUTLINE OF THE THEORY 

A. RSJ + Johnson noise 

The simplest model to describe thermal phase fluctua- 
tions in Josephson systems is the usual RSJ model with 
Johnson noise 



fiC d if h dip 
2e dt 2 + 2eRdi 



(Mt)J T (t')) 



J c sin ip = J + Jt, 



2kT 



5(t-t'). 



(1) 



(2) 



Here ip is the Josephson phase difference, C, R(T), J C {T) 
are junction capacitance, resistance, and critical current, 
J is the bias current, Jr(t) is a random current noise. 
This equation can be written in dimcnsionlcss form as 



d 2 ip dip 



+ — + smip = j + j T , 



dr 2 dr 



(3) 



As mentioned above, for the description of our experi- 
ments we need a frequency dependent damping which is 
achieved by the coupling to an external circuit as shown 
in Fig. |21 This model was already considered in detail by 
Kautz and Martinis* and is described by the following 
set of equations 



dr 2 



dip 



dp 



+ sin ip = j + j T + j Ts , (7) 



pdv = /dip 
dr \ dr 



JTs 

a 



where 



(8) 



(9) 



Here the parameters a = R/R s and p = aC/C s are intro- 
duced, v(t) is the voltage across the external capacitance. 

From the geometry of our samples we assume that C s 
is the capacitance between the top and bottom contacts 
of the mesa and consequently C s ^ C . The external 
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Figure 3: Typical current-voltage curves at low temperatures. 
The zero temperature critical current is suppressed by prema- 
ture switching 
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Figure 4: Typical current-voltage curves at intermediate tem- 
peratures with frequency-dependent damping (circles) and in 
the simple RSJ model (triangles). 



resistance R s is determined by a contact resistance and 
the resistance of the leads. For small mesas the single 
junction resistance R is very large, up to 10-100 kf2 and 
fl s «ii can be assumed. Therefore for small junctions 
the coupling to the external circuit should be included 
in the model calculations. In the classical fluctuation 
regime discussed in our paper it leads to the existence of 
a phase-diffusion branch instead of a zero-voltage super- 
conducting branch. 

Using this model as basis for numerical simulations we 
discuss in the following the typical behaviour of current- 
voltage curves in the different temperature regimes. In 
Section IV the model will be refined, in order to describe 
the experimental results in detail. 



C. Premature switching at low temperatures 

(T < T*). 

At low temperatures a typical I-V curve as shown in 
Fig-El contains a superconducting branch (S-state) with 
zero voltage and a hysteretic jump to the resistive state 
(R-state). Due to fluctuations the critical (or, better 
to say switching) current J s is random: the first strong 
enough thermal kick leads to a premature transition to 
the R-state, the probability of a back-transition to the 
S-state is very small. The average switching current 
(J S )(T) depends on the sweeping velocity I of the cur- 
rent. It can be calculated numerically and analytically. 
We found that the results of our numerical simulations 
are in good agreement with the well known analytical 
result (J C )(T) at low temperatures^ 



<J S )(T) = J C (T) 



2 27r/ 



2/3' 



(10) 



Typically the logarithm is of the order of 10. 



D. Phase diffusion branch at intermediate 
temperatures (T ~ T*) 

At intermediate temperatures the S-branch shows a 
bending with a finite voltage close to the switching cur- 
rent (see Fig.0J). As mentioned above this behaviour 
can be explained by phase diffusion within a model with 
frequency-dependent damping. Let us discuss briefly the 
properties of this model. 

In the limit of zero frequency dissipation is determined 
by the junction resistance R alone, but at frequencies of 
the order of the plasma frequency oj p due to the small 
impedance 1 /oj p C s of the external capacitance the dissi- 
pation is determined by the external resistance R s <C R 
in parallel with R. Therefore for small junctions with 
R s <C R dissipation at the plasma frequency is signif- 
icantly larger then at zero frequency. For large-/? hys- 
teretic junctions it leads to the coexistence of a finite- 
voltage phase-diffusion S-branch and resistive branch. 
Thermal current-kicks stimulate jumps of the Josephson 
phase difference from one potential minimum to the next 
and therefore diffusion. On the other hand, the high 
dissipation at the plasma frequency - the characteristic 
frequency of these jumps, prevents from a fast transition 
to the resistive state with a running phase. 

In Fig. 01 the phase diffusion branch in the model with 
frequency-dependent damping is shown (circles) together 
with the usual RSJ curve at similar parameters (trian- 
gles). 

The phase diffusion branch can be found numerically, 
and has been observed in our experiments, only at inter- 
mediate temperatures. At low temperatures the voltage 
in the S-state is exponentially small. At high tempera- 
tures the usual phase diffusion takes place with a single- 
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Figure 5: Typical current-voltage curves at high tempera- 
tures. The switching current J a is determined as the current J 
for which the contribution of an average supercurrent (dashed 
line) is maximal 



valued voltage-current curve. 



E. Switching current at high temperatures 

(T > T*). 

At high temperatures the hysteresis disappears and 
only one single-valued I-V curve with a finite voltage at 
all currents exists (Fig.|SJ, but the average supercurrent, 
which is large at low voltages leads to a pronounced non- 
linearity of I-V curve. We can extract approximately the 
contribution of the supercurrent to the I-V curve by sub- 
tracting the quasi-particle current, which in this model is 
a linear function of voltage. The current J at which the 
difference J — V/R reaches its maximum (see Fig.[5J) de- 
fines the switching current J s . Note that this procedure 
to determine a switching current can also be applied at 
low temperatures and will be generally used later when 
experimental and theoretical results are presented and 
compared. 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL SET-UP 



To measure the properties of ultrasmall intrinsic 
Josephson junctions we used a mesa geometry as sketched 
in Fig.© Bi 2 Sr 2 CaCu 2 8+5 (BSCCO) single crystals 
and Tl 2 Ba 2 CaCu20 8+ 5 (TBCCO) thin films were used 
to prepare the samples^. Fig. shows an SEM pic- 
ture of a TBCCO sample. With our technique we were 
able to structure samples with lateral dimensions as small 
as 0.2 x 0.2/jto 2 . The height of the mesas was between 
60 A and 1500 A, which is equivalent to 4 - 100 intrinsic 
Josephson junctions. 

Several small mesas with lateral dimensions of a few 
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Figure 6: Sketch of the sample geometry 



Small mesa 




Large mesas 

Figure 7: SEM picture of a TBCCO sample. 



\xm or lower were arranged on one chip. On the same chip 
also some larger mesas with lateral dimension of several 
/im were prepared. They serve as electrical contacts to 
the bulk superconductor. All mesas were contacted with 
Au leads. Current flows through one large mesa into the 
bulk superconductor and is extracted through the small 
mesa or vice versa. For voltage measurements one of the 
other large mesas is used. Due to the dimensions of the 
small mesa it was not possible to attach two isolated Au 
leads. Thus transport measurements were performed in 
a three-terminal configuration. 

The Au electrode connected to the mesa leads to a sup- 
pression of the critical current and critical temperature 
of the uppermost intrinsic Josephson junctioniiiS. This 
degeneration might be explained by the assumption that 
the first superconducting layer is in close proximity with 
the normal Au electrode. 

One sample was prepared in a step like geometry as 
proposed by Wang et ali&. Here we are able to record the 
/-^-characteristics in a true four-point geometry (sample 
s7n4 jjl Tabled). 

The bias current was provided by a battery powered 
current source. The 7-V~-characteristics of the samples 
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Table I: Sample parameters. J s is the switching current. T* xp 
is defined by the kink in the experimental J s vs. T plots (see 
Fig.0 




were recorded by applying dc currents and recording the 
voltages across mesas by digital voltmeters. 

Measurements were done either in a standard helium 
dewar or in a magnet cryostat equipped with a 5 T 
Helmholtz solenoid. Temperatures could be varied be- 
tween 4.2 K and room temperature. The temperature 
was measured with a platinum resistor for higher tem- 
peratures and a cernox resistor for low temperatures. 

To reduce external noise low pass filters were used. In 
the helium dewar both cold filters and room temperature 
filters were used. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Experimental voltage-current curves and J S {T) 

We systematically measured voltage-current curves of 
mesas with different zero-temperature critical currents 
J c (0) at different temperatures (only a part of the sam- 
ple parameters is presented in Table. I). Transition from 
high-hysteretic curves at low temperatures to phase- 
diffusion branches at crossover temperatures, and to a 
single resistive curve at high temperatures was clearly ob- 
served. Examples of experimental current-voltage curves 
at low (squares) and intermediate (filled circles) temper- 
atures are shown in Fig-El (here the contact resistance is 
not subtracted). The well-known multibranch structure 
is observed. 

The measurements of critical current were done mostly 
for the first (superconducting) branch, such that all junc- 
tions are in the S-state at low current and only one junc- 
tion with the lowest critical current switches to the Re- 
state. For this branch the average switching current J a 
was determined. Its temperature dependence is shown in 
Fig-Elfor five samples with different values of J c (0) (open 
symbols) together with the results of the theoretical cal- 
culation (filled symbols, discussed in the next section). 
This picture is the main result of our work. 

The main feature of all curves is a kink (or even min- 
imum) in J S (T) at some temperature T* . This kink is 



Figure 8: Typical experimental current-voltage curves at low 
(squares) and intermediate temperatures (circles). 
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Figure 9: Js(T) for samples with different J c (0). Experiment 
(open symbols, the lines are guides to the eye), and theory 
(filled symbols). 



absent in the J S (T) curves of mesas with large values of 
J c (0) and is shifted to lower temperatures with decreas- 
ing J c (0). T* xp corresponds approximately to the con- 
dition Ej(T*) — kT*, which strongly suggests thermal 
fluctuations as an origin of the observed features. 



B. Theoretical voltage-current curves and Js(T) 

To describe the experimental result for the switching 
current J s (T) in detail it is necessary to take into account 
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Figure 10: Quasiparticle conductivity of an intrinsic Joseph- 
son junction as function of T shown for incoherent tunneling 
and unitary intra- layer scattering (solid line), experimental 
points (crosses), and the fit g(T) used in further calculations 
(circles). 
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Figure 11: Phase-diffusion at low currents. Experimental 
points at different temperatures and fit to the analytical ex- 
pression itTTl . The imped ance Z\ (shown in the insert) was 
used as a fitting parameter. 



the temperature dependence of the parameters f3 and 7. 
These parameters are known functions of: 

(i) the critical current J C (T) without fluctuations, 
which is determined (qualitatively also in HTSC) by the 
Ambegaokar-Baratoff relation. 

(ii) the quasiparticle resistance R(T), which we de- 
termined experimentally from current-voltage curves at 
small voltages (crosses in Fig. llOfl . We also calculated 
R(T) at V — » from microscopic theory for different 
models of interlayer tunneling (coherent and incoherent) 
and intra-layer scattering (Born and unitary). The best 
fit to the experimental results was obtained for incoher- 
ent tunneling in the unitary limit of intra-layer scattering 
with an elastic scattering frequency v = 1 meV. For the 
numerical simulations the function g(T) = R(0)/R(T) is 
approximated by the fit (circles in Fig.llO[). 

The results of the theoretical calculations are presented 
in Fig.[5]for three different zero-temperature critical cur- 
rents J c (0), the zero-temperature McCumber parameter 
(3(0) is in the range of 10 3 to 10 4 . For each tempera- 
ture we calculated J C (T), (3(T), j(T), and then produced 
current-voltage curves by numerical simulations starting 
from zero external current J. From these results the cur- 
rent J was determined, at which the average supercur- 
rent J SU p = J — V/ R reaches its maximum. This current 
was defined previously as switching current J s . The ex- 
perimental switching currents J s shown in Fig.[5]are de- 
fined in the same way. We find good agreement between 
experimental and theoretical curves. Also the current- 
voltage curves calculated at different temperatures show 
the same behaviour as the experimental curves. 

From these results we conclude that the non-monotonic 
temperature dependence of J s (T) observed in experiment 
can well be explained by the model of large-/? Josephson 
junction with thermal noise and coupling to an exter- 
nal impedance, if the proper temperature dependence of 



junction resistance is taken into account. In particular 
the strong suppression of J s with temperature at low tem- 
peratures, T < T* , is due to the premature switching ef- 
fect in highly-hysteretic junctions. At high temperatures, 
T > T*, the McCumber parameter (3 is smaller, the hys- 
teresis disappears and the contribution of the average 
supercurrent to the current at finite voltages becomes 
larger, which leads to the plateau in the observed J S {T) 
dependence. Further discussions of the origins of non- 
monotonic temperature dependence of J s can be found 
in the papers of Iansiti et a\& and Kautz and Martinis 7 . 



C. Phase-diffusion at low currents 

For the junctions with the smallest critical cur- 
rent J c (0) ~ 1 fiA and highest quasiparticle resistance 
R^> R s we investigated the low-current part of the 
current- voltage curves more precisely and found clear ev- 
idence of the phase-diffusion behaviour induced by the 
external circuit. 

In the case R s <C R and for Ej < kT the phase- 
diffusion branch can be described by the analytical 



expressio: 



.15.20 



J{V) 



4eTJ, 



Z X V 



h V 2 + {2eZ l T/hf 



(11) 



and the small finite resistance in the S-state at very low 
currents and at any temperature is given bySS 



R( J ->0) 



Wr)-i' 



(12) 



where Z\ ~ R s is the impedance at low frequencies (of 
the order of the plasma frequency u> p ) and Iq(x) is a 
Besscl function. 
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Figure 12: Js(T) at different perpendicular (upper figure) and 
parallel (lower figure) magnetic fields: experiment. 



dure described in Ref. [^J The junction is biased at a 
fixed current and a fixed magnetic field is applied. Then 
the voltage change is monitored while the angle cp be- 
tween the external field and the Cu02-layers was var- 
ied with a precision rotation stage. When measured at 
high enough temperatures the field-induced changes of 
the /-^-characteristic were reversible and could be used 
for alignment. Thus, the field orientation relative to the 
CuCVlayers could be adjusted to an accuracy of 0.01°. 

The main and most general experimental result is pre- 
sented in Fig. 1121 In both perpendicular and parallel 
magnetic fields J S (T) is not changed for T > T*, while 
for T < T* one finds a field dependent suppression of 
J S {T). 

To understand this behaviour qualitatively, we explore 
the well-known sin-Gordon equation, which describes 
long Josephson junctions (see Ref. l2^l23|) 

d 2 ip d 2 ip dtp 
/3-^--j-^- + -^:+sin(<^+ J B||X+(5^(a;,T)) = j+j T (x,r). 

(13) 

Here the coordinate x is in units of the Josephson length 
A j, B\\ is the component of the external magnetic field 
parallel to the layers, and Stp(x, r) is a random phase dis- 
tribution produced by pancake-vortices in perpendicular 
magnetic field24iSiS&. We include the phase variation due 
to the external magnetic field in the sin-function instead 
of using the more familiar procedure, where the magnetic 
field is taken into account by boundary conditions (both 
methods are equivalent). The random functions jx(x,T) 
and 5ip(x, r) describing noise are assumed to be Gaussian 
and 5-correlated in space 

(jT(x,T)j T (x,T)} = 2-f5(r - t')5(x - a;'), (14) 



The experimental points presented in Fig. 1111 can be 
fitted by Eq. The best fitting result was achieved 

for a temperature dependent impedance Z\ as shown in 
the insert. The origin of this temperature dependence is 
still unclear. It may be a contribution of the temperature 
dependent quasiparticle conductivity of other junctions 
in the mesa: when the temperature increases this con- 
ductivity also increases and the impedance decreases. 



D. Js(T) in external magnetic field 

Finally, we present here some results of measurements 
in an external magnetic field. Magnetic fields were ap- 
plied by a 5 T superconducting split coil magnet. To 
measure J S (B) for parallel fields, the Cu02-layers must 
be aligned with high precision parallel to the magnetic 
field. If the magnetic field has a component perpendic- 
ular to the layers pancake vortices enter the Josephson 
junctions and suppress the critical current. The orien- 
tation of samples in a magnetic field follows the proce- 



(5ip(x, t)5<p(x 7 t')) oc B z S(x - x'). (15) 

In the last expression we accept a linear magnetic field 
dependence of the correlation function and neglect the 
time-dependence of phase shifts assuming that the pan- 
cakes are frozen or move slowly enough. In a more real- 
istic approach the space and time dependence of correla- 
tion functions of the pancake disorder depend on a state 
of the vortex matter, is it liquid, glass-like, or solid. 

Results of our numerical simulations with magnetic 
field are presented in Fig.^| for several values of per- 
pendicular and parallel fields. The main result of our 
calculations (as well as similar experimental results) is 
the sensitivity to a magnetic field at low temperatures 
and the absence of considerable influence at high tem- 
peratures. 

Qualitatively this result can be understood in the fol- 
lowing way. The main effect of the magnetic field is to 
produce phase shifts, regular for the parallel field and 
random for the perpendicular field. At low tempera- 
tures, when J s is determined by the premature switching 
mechanism, the reduction of J c in the magnetic field is 
important because the energy barrier becomes smaller. 
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Figure 13: Js(T) at different perpendicular (circles), parallel 
(squares), and oblique (triangles) magnetic fields: theory. $ 
is the magnetic flux of the parallel field through a Josephson 
junction in units of the flux quantum $o- 



But in the high-temperature regime where phase diffu- 
sion is strong, additional phase shifts 5<p(x, r) produced 
by the magnetic field cannot compete with the existing 
large thermal phase fluctuations. 



In this paper we discuss new experiments on the in- 
fluence of thermal noise on current-voltage curves and 
the switching current of small area intrinsic Josephson 
junctions. Two main effects of thermal fluctuations: pre- 
mature switching and phase-diffusion, are observed. The 
non-monotonic temperature dependence of the switch- 
ing current is explained theoretically with a help of an 
extended RSJ model (including thermal noise and a fre- 
quency dependent external impedance) taking into ac- 
count the proper temperature dependence of parameters. 
The different behaviour of the switching current at low 
and high temperatures in the presence of magnetic fields 
can be explained in a similar way. 
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